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摘 要 : Sentinel- A 合成 孔径 雷达 数据 不 受 云雾 等 天 气 条 件 的 影响 并 具有 丰富 的 纹理 信息 ,为 提 


取 城 市 建设 用 地 信息 提供 了 一 种 新 的 数据 源 。 本 文 发 展 了 


种 基于 Sentinel-1A 合成 孔径 雷达 数 


据 和 全 卷 积 网 络 的 城市 建设 用 地 监测 方法 。 该 方法 的 优势 主要 在 于 可 以 有 效 复合 不 同 极 化 方式 下 
的 Sentinel-1A 合成 孔径 雷达 数据 和 综合 集成 多 尺度 特征 。 在 甘肃 省 张掖 市 甘 州 区 的 应 用 表明 ,该 
方法 的 提取 结果 总 体 精 度 为 92.50% ,Kappa 系数 为 0.85。 与 现 有 方法 KTH-Pavia 城市 提取 器 ” 相 
比 ,Kappa 系数 提高 了 37.10% ,总 体 精 度 提高 了 11.50% 。 因 此 ,该 方法 具有 良好 的 应 用 潜力 。 
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城市 建设 用 地 是 指 包括 居住 用 地 、 商 业 服 务 业 
设施 用 地 、 工 业 用 地 .道路 交通 设施 用 地 和 公用 设施 
用 地 在 内 的 城市 土地 ”。 伴 随 着 经 济 的 高 速 发 展 
与 城市 人 口 的 持续 增长 ,中 国正 在 经 历 大 规模 的 城 
市 建设 用 地 扩展 过 程 *”“ 。 该 过 程 一 方面 促进 了 社 
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刘 纪 远 等 2” 以 Landsat 数据 为 主要 数据 源 ,获取 了 
中 国 20 世纪 80 年 代 ~ 2010 年 的 城市 建设 用 地 信 
Ho HE 4°! 复合 MODIS 数据 以 及 DMSP-OLS、 
NPP-VIIRS 夜间 灯光 数据 获取 了 中 国 1992 一 2012 
年 城市 建设 用 地 信息 。 但 光学 遥感 易 受 云 . 雾 等 天 


会 发 展 并 提高 了 城市 居民 的 人 类 福 社 , 另 一 方面 也 
带 来 了 水 资源 短缺 .空气 污染 和 城市 热岛 效应 等 一 
系列 问题 “ … 。 及 时 准确 地 提取 城市 建设 用 地 信 
息 是 理解 并 评价 城市 扩展 过 程 对 生态 环境 影响 的 基 
础 ,对 于 促进 区 域 可 持续 发 展 具有 重要 意义 ”"]。 
光学 遥感 数据 是 当前 获取 城市 建设 用 地 信息 的 
主要 数据 源 。 包 括 Landsat 数据 .中 巴 资源 卫星 数 
据 和 高 分 系列 卫星 数据 等 中 高 分 辩 率 遥感 数据 ,以 
及 MODIS ( Moderate Resolution Imaging Spectroradi- 
ometer ) 数据 、DMSP-OLS Defense Meteorological Sat- 
ellite Program’s Operational Linescan System ) 夜间 灯 
光 数 据 和 NPP-VIIRS( National Polar-orbiting Partner- 
ship’s Visible Infrared Imaging Radiometer Suite ) 夜间 
灯光 数据 等 中 低 分 状 率 遥感 数据 在 内 的 光学 遥感 数 
据 已 被 广泛 用 于 城市 建设 用 地 的 提取 "1。 例 如 ， 
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气 条 件 影响 ,同时 其 所 包含 的 城市 地 物 纹理 信息 不 
Æ ,在 用 于 获取 城市 建设 用 地 信息 时 仍 受到 限制 。 
合成 孔径 雷达 (synthetic aperture radar，SAR ) 
数据 为 提取 城市 建设 用 地 信息 提供 了 一 种 新 数据 
Wo SAR 是 一 种 高 分 辩 率 微波 成 像 雷 达 , 主要 原理 
是 通过 小 孔径 雷达 天 线 沿 直线 不 断 移动 ,在 不 同位 
置 上 接收 同一 地 物 的 回 波 信 号 并 通过 数据 处 理 合 
成 ,从 而 等 效 为 一 个 大 孔径 天 线 ,实现 高 分 辨 率 成 
RUH, SAR 有 着 全 天 候 、 全 天 时 工作 能 力 ,雷达 天 
线 所 接收 的 雷达 波 为 后 向 散射 信号 。 在 城市 中 , 平 
顶 建 筑 .街道 等 对 雷达 波 产生 镜面 反射 , 相 较 周 围 地 
物 呈 暗色 调 ,建筑 物 的 立体 结构 则 会 产生 “* 角 反射 
器 "效应 增强 回 波 信号 ,从 而 在 影像 上 表现 出 丰富 
的 纹理 特征 。 此 外 ,数据 的 多 极 化 类 型 可 以 反映 
散射 体 结构 。 因 此 ,雷达 数据 适合 用 于 城市 建设 用 
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地 的 提取 “-” 。 已 有 相关 研究 所 采用 的 SAR 数据 
主要 包括 Radarsat-208 -2 、ENVISATI2 -° ERS - 
1/ -2!"-73) 和 TeraSAR-X2 >], {B ERS - 1/-2, 
ENVISAT 等 早期 卫星 已 停止 工作 ,TeraSAR-X、Ra- 
darsat-2 等 获取 的 数据 价格 昂贵 ,限制 了 SAR 数据 
的 进一步 应 用 ”1 。 

Sentinel-1A 数据 因 具 有 免费 .便捷 获取 等 优势 ， 
已 成 为 一 种 被 广泛 应 用 的 SAR 数据 源 '” o, 2014 年 
4 月 3 日 发 射 的 Sentinel-1A 卫星 搭载 了 C 波段 合成 
孔径 雷达 , 重 访 周 期 为 12 dL28-29] 。 Sentinel-1A 数据 
包含 多 种 极 化 类 型 ,能 以 较 大 幅 宽 实现 全 球 陆地 履 
盖 , 兼 有 重 返 周 期 短 和 免费 获取 等 突出 优势 ,在 城市 
建设 用 地 提取 研究 中 被 广泛 应 用 中。 例如 ,BAN 
等 所 以 Sentinel-1A 为 主要 数据 源 ,在 全 球 范围 内 提 
取 了 十 个 代表 性 城市 的 建设 用 地 信息 。 但 基于 Sen- 
tinel-1A 数据 提取 城市 建设 用 地 的 现 有 方法 还 存在 
AR a AS A) Aan, GAMBA 40-71 采用 的 
KTH-Pavia 城市 建设 用 地 提取 方法 主要 通过 复合 空 
间 自 相关 指数 与 灰 度 共生 算 阵 提取 城市 建设 用 地 信 
息 ,该 方法 未 能 充分 融合 Sentinel-1A 数据 的 不 同 极 
化 类 型 和 不 同 尺度 的 纹理 信息 , 阅 值 的 选取 也 比较 
主观 ,提取 效果 有 较 大 提升 空间 。 

全 卷 积 网 络 提 供 了 提取 城市 建设 用 地 信息 的 新 
途径 。 全 卷 积 网 络 是 LONG 4] 于 2015 年 在 卷 积 
神经 网 络 基础 上 发 展 出 的 一 种 采取 “ 像 元 到 像 元 ” 
图 像 识别 模式 的 深度 学 习 网 络 结构 。 全 卷 积 网 络 主 
要 由 卷 积 层 和 激活 层 组 成 ,继承 了 卷 积 神经 网 络 的 
基本 结构 。 其 中 , 卷 积 层 能 够 识别 图 像 局 部 特征 , 激 
TREY LIA TATE ERE OP, HL BAS 
积 层 和 激活 层 交 替 传 递 信息 ,全 卷 积 网 络 能 有 效 复 
合 多 源 信 息 ,综合 集成 多 尺度 特征 ,从 而 对 地 物 进行 
准确 分 类 “]。 近 年 来 ,已 有 研究 人 员 利 用 全 卷 积 
网 络 开 展 了 基于 光学 遥感 影像 的 城市 建设 用 地 提取 
研究 , 且 取 得 了 良好 的 效果 。 例 如 ,FU EY 利用 全 
卷 积 网 络 和 GF-2 数据 提取 了 北京 东北 部 城市 建设 
用 地 ; ZHANG 等 “利用 全 卷 积 网 络 和 WorldView 
高 分 影像 数据 获取 了 北京 市 海淀 区 与 西城 区 城市 建 
设 用 地 信息 。 但 利用 全 卷 积 网 络 和 Sentinel-1A 数 
据 提 取 城 市 建设 用 地 信息 的 研究 方法 仍 待 探索 。 

本 文 则 在 发 展 一 种 利用 全 卷 积 网 络 的 城市 建设 
用 地 监测 方法 ,基于 Sentinel-1A 合成 孔径 雷达 数据 
实现 城市 建设 用 地 的 有 效 提取 。 该 方法 不 仅 能 充分 
发 挥 Sentinel-1A 合成 孔径 雷达 数据 的 优势 ,同时 具 


备 良好 的 应 用 潜力 。 
1 研究 区 概况 


研究 区 位 于 中 国 甘肃 省 张掖 市 甘 州 区 ,地 理 位 
置 为 100°22’15” ~ 100°32'20"E ,38°54'31” ~ 38°58’ 
54"N( 图 1)。 研 究 区 地 物 中 的 裸 地 与 城市 建设 用 地 
光谱 特征 相近 , 易 导 致 混 分 问题 。 而 张掖 市 是 中 
绿洲 城市 土地 利用 /覆盖 变化 研究 的 热点 城市 “。 
在 “西部 大 开发 ”等 政策 的 推动 下 , 随 着 经 济 和 人 口 
规模 的 高 速 增长 ,该 区 域 近年 来 经 历 了 快速 的 城市 
建设 用 地 扩张 过 程 ”” 。 综 上 ,研究 区 适 于 开展 基 
于 Sentinel-1A 合成 孔径 雷达 数据 的 城市 建设 用 地 
提取 方法 研究 。 


2 数据 


所 使 用 的 Sentinel-1A 合成 孔径 雷达 数据 来 源 
于 欧洲 航空 局 哥 白 尼 全 球 对 地 观测 项 目 公 开发 布 的 
Sentinel-1 A 数据 集 ( https://scihub. copernicus. eu/ )。 
Sentinel-1A 合成 孔径 雷达 包括 4 种 工作 模式 ,分 别 
是 条 带 模 式 (Strip Map Mode,SM) .干涉 宽 幅 模式 
(Interferometric Wide Swath ,IW ) 、 超 宽 幅 模式 ( Extra 
Wide Swath EW) \ 波 模式 (Wave Mode,WV) ”。 根 
据 数 据 的 可 获取 性 与 空间 分 辨 率 , 选 择 使 用 IW 模 
式 数据 ,该 模式 为 卫星 默认 的 对 地 观测 模式 , 幅 宽 
250 km ,空间 分 辨 率 为 5 m x20 m, 极 化 方式 包括 
VH( 交 叉 极 化 ) 与 VV( 垂 直 同 极 化 )2 BROT, BHR 
获取 时 间 为 2018 年 2 月 5 日 ,空间 分 辩 率 重 采 样 为 
10 m x 10 m。 参 考 BAN 等 的 研究 ,经 过 辐射 校 
IE .斑点 噪声 抑制 和 几何 校正 等 预 处 理 后 ,将 DN 值 
标准 化 到 0 ~ 255 之 间 。 合 并 2 种 极 化 方式 为 双 波 
段 影 像 数据 ,VH 极 化 数据 作为 第 一 波段 , VV 极 化 
数据 作为 第 二 波段 。 

用 于 获取 训练 样本 和 精度 评价 的 高 分 辨 率 遥 
感 影像 数据 来 源 于 Coogle Earth 平台 提供 的 Digital 
Globe 3 4% (http://www. digitalglobe. com/) ,空间 
分 辩 率 约 为 7.4 吓 ,影像 拍摄 时 间 为 2018 年 3 月 7 
日 。 基 于 高 分 影像 数据 进行 人 工 目 视 解 译 分 类 出 
城市 与 非 城市 训练 样本 。 为 确保 精度 ,于 2018 年 
7 月 1~2 日 对 研究 区 开展 了 实地 调查 ,结合 调查 
结果 对 目 视 解 译 结果 进行 修正 。 最 终 获得 的 地 面 
真 值 用 于 全 卷 积 网 络 的 训练 和 提取 结果 的 精度 
评价 。 
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Sentinel- 1A 合 成 孔径 雷达 数据 (VH 极 化 ) 


到 1 研究 区 


Fig.1 The study area 


(a) 整体 技术 路 线 
'Sentinel-1A 合 
成 孔径 雷达 数据 
b) 构 wensenmed _ 


eet 层 激 活 层 


图 2 技术 路 线 
Fig.2 Flowchart 
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3 方法 


3.1 构建 全 卷 积 网 络 

参考 LONG 28) Al SZEGEDY 22°! 的 研究 , 构 
建 了 一 个 全 卷 积 网 络 用 以 提取 城市 建设 用 地 信息 
(图 2)。 该 网 络 的 基本 思路 是 利用 输入 层 复合 不 同 
极 化 方式 的 Sentinel-1A 数据 ,而 后 利用 基于 Incep- 
tion 结构 "的 卷 积 层 .激活 层 和 连接 层 获取 不 同 尺 
度 的 纹理 特征 ,最 后 通过 分 类 器 提取 出 城市 建设 用 
地 信息 (图 2a)。 

为 了 确定 卷 积 核 的 大 小 和 卷 积 层 的 层 数 ,利用 
验证 集 对 基于 不 同 大 小 卷 积 核 和 不 同 层 数 卷 积 层 构 
建 的 网 络 进行 了 测试 ,根据 分 类 精度 和 运行 效率 ,最 
终 选 择 使 用 由 五 个 Inception 结构 组 成 的 深度 全 卷 
积 网 络 结构 (图 2b)。 模 型 的 输入 层 接收 包括 VH 
和 VV2 种 极 化 方式 的 Sentinel-1A 数据 。 卷 积 层 的 
3 种 卷 积 核 大 小 分 别 为 1 x1、3 x3 和 5 x5, 用 于 获 
取 不 同 尺度 的 特征 信息 。 卷 积 层 与 输入 层 的 关系 可 
用 下 式 表示 。 


Conv; = so we x Input’, ;) + Bias” (1) 
RP : Conny) RAN | MERRY r 种 卷 积 核 中 第 
i 个 波段 第 7 个 像 元 的 DN 值 ;Inputi ,表示 第 1 层 输入 
数据 中 第 了 个 波段 第 7 个 像 元 的 DN 值 ;Wi 表示 该 
卷 积 层 的 权 值 ; Bias;” 表示 偏 置 项 ; No 表示 输入 数 
据 的 波段 数 ;K, 表示 第 了 种 卷 积 核 的 大 小 。 连 结 层 
用 于 融合 上 一 个 卷 积 层 中 3 种 卷 积 核 所 提取 的 信 
息 ,连结 层 的 输出 结果 作为 下 一 个 卷 积 层 的 输入 数 
据 。 连 结 层 与 上 一 层 卷 积 层 的 关系 可 用 下 式 表 示 。 


f( Conv; ) i<m 


Cone; ; = f Con; nj) m<i<me+n (2) 


及 Conv? n -nj ) 


式 中 ;Conci ,表示 第 /个 连结 层 中 第 i 个 波段 第 j 个 
像 元 的 DN 值 ;m,n,k 分 别 代 表 该 层 三 种 卷 积 核 的 
数量 ;f(x) 表 示 激 活 函 数 ,所 使 用 的 激活 函数 是 标准 
的 修正 线性 单元 ( Rectified Linear Unit,ReLu) ,可 表 
示 为 : 


m+n<i<m+nt+tk 


f(x) =max(0,x) (3) 
设置 多 个 Inception 结构 的 目的 是 增 大 有 效 感 
受 野 ,提高 全 卷 积 网 络 对 多 尺度 特征 信息 的 识别 能 


力 。 最 后 两 个 卷 积 核 获 得 各 像 元 被 分 类 为 城市 和 非 
城市 的 概率 ,可 用 下 式 表示 。 


1, Output” > Output?" ae 
Urban, = | (4) 
0, Otherwise 
SUP: Urban, 表示 输出 层 中 第 7 个 像 元 是 否 为 城市 
建设 用 地 ,1 为 城市 建设 用 地 ,0 为 非 城市 建设 用 地 ; 
Output" 表示 第 j 个 像 元 为 城市 建设 用 地 的 概率 ; 
Output?" 表示 第 j 个 像 元 为 非 城市 建设 用 地 的 
概率 ,可 通过 下 式 计算 。 


k 
Output" = f( > We Cy x Cone; ; ) + 


Bias urban ) (5 ) 


m+n+k 


Output" = f( > 


( We x Cone; ,) i 


i 
Bias’ ) (6 ) 


式 中 : wri 和 Ween Re A LAL Bias” 和 
Bias "表示 偏 置 项 。 上 述 全 卷 积 网 络 的 构建 依 
托 Convolutional Architecture for Fast Feature Embed- 
ding( Caffe) REED EASY 。 
3.2 训练 全 卷 积 网 络 

训练 全 卷 积 网 络 的 目的 是 求 算 所 构建 的 全 卷 积 
网 络 中 权 值 的 全 局 最 优 解 。 基 本 原理 是 以 损失 函数 
最 小 为 求解 目标 ,利用 优化 算法 迭代 求解 最 优 权 值 。 
步 又 主要 包括 选择 训练 样本 确定 损失 函数 和 设置 
优化 算法 `. 友 代 参数 三 个 基本 步骤 。 首 先 ,以 Google 
Earth 高 分 影像 数据 为 基础 ,通过 目 视 解 译 提取 出 城 
市 与 非 城 市 训练 样本 图 像 , 并 结合 野外 调查 进行 修 
正 。 将 整个 训练 样本 图 像 裁 前 为 32 x 32 像素 大 小 
的 图 像 块 ,共计 1 534 块 。 按照 8:2 比例 随机 划分 
训练 集 与 验证 集 , 获 得 训练 集 1 226 块 ,验证 集 308 
块 。 

其 次 ,参考 LONG 等 ”的 研究 ,在 全 卷 积 网 络 
的 输出 层 后 加 入 用 于 计算 损失 函数 的 SoftmaxWith- 
Loss 层 。 损 失 函 数 可 用 下 式 表示 。 


I(y,z) =log Ye’ -z, (7) 


式 中 :z 表示 SoftmaxWithLoss 层 的 输入 ;1(y,z) 表 示 
z 的 损失 函数 ;sz 表示 第 7 个 类 别 的 预测 结果 ;y 表示 
训练 样本 中 z 的 类 别 。 

最 后 ,参考 LONG 等 中 的 研究 ,确定 了 优化 算 
法 与 迭代 参数 。 优 化 算法 使 用 小 批 次 随机 梯度 下 降 
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法 (mini-batch Stochastic Gradient Descent ) , 即 在 每 
次 迭代 中 使 用 不 同 的 小 批 次 训练 数据 计算 后 向 传播 
梯度 ,并 在 近 代 过 程 中 不 断 减 小 参数 修正 值 直 到 将 
损失 降 至 最 低 。 超 参数 是 根据 经 验 值 和 模型 实际 训 
练 效果 获取 ,初始 学 习 率 (base_T) 设 为 8 x 107° HE 
次 大 小 (batch size) 设 为 8 ,学 习 率 下 降 方式 (Jr_poli- 
cy) 为 “step”, 迭 代步 长 (step size) 设 为 5 000 ,衰减 
量 (Gamma) 设 为 0.99 ,动量 (Momentum ) 设 为 0.99。 
利用 训练 后 的 全 卷 积 网 络 ,基于 Sentinel-1A 数据 提 
取出 了 研究 区 2018 年 城市 建设 用 地 信息 。 


4 结果 


为 检验 本 方法 的 提取 效果 ,参考 JANSSEN 等 ”| 
的 研究 ,采用 分 层 随 机 抽样 的 样本 点 检验 方法 ,在 整 
个 研究 区 分 别 选 取 了 100 个 城市 样本 点 和 100 个 非 
城市 样本 点 ,结合 同期 高 分 状 率 遥感 影像 ,通过 构建 
混淆 矩阵 并 计算 Kappa 系数 ,总体 精度 .制图 精度 和 
用 户 精 度 等 主要 评价 指标 对 提取 结果 进行 精度 评价 
( 表 1) 。 结 果 表 明 ,基于 全 卷 积 网 络 提取 的 城市 建设 
用 地 Kappa 系数 达到 0. 85 ,总 体 精 度 .制图 精度 与 用 
户 精 度 均 在 90% 以 上 ,说 明 提 取 结 果 可 以 准确 可靠 


表 1 基于 全 卷 积 网 络 的 城市 建设 用 地 提取 结果 的 精度 评价 
Tab.1 Accuracy assessment of urban construction land 
extraction result based on fully convolutional network 

参考 图 像 
城市 建设 非 城市 建设 ”总 和 
用 地 ( 像 元 ) 用 地 ( 像 元 )( 像 元 ) 


城市 建设 用 地 ( 像 元 ) ”96 4 100 
分 类 图 像 ” 非 城市 建设 用 地 ( 像 元 ) 11 89 100 
总 和 ( 像 元 ) 107 93 200 


总 体 精 度 =92.50% Kappa =0. 85 
制图 精度 ( 漏 分 误差 ) 用 户 精 度 ( 错 分 误差 ) 
89.72% (10. 28% ) 96.00% (4.00% ) 
95.70% (4. 30% ) 89.00% (11.00% ) 


城市 建设 用 地 
非 城市 建设 用 地 


地 反映 研究 区 城市 建设 用 地 的 空间 格局 。 

提取 结果 表明 ,张掖 市 甘 州 区 2018 年 建设 用 地 
面积 为 4 006 hm , 约 占 区 域 总 面积 的 20% (图 3)。 
城市 建设 用 地 斑 块 的 面积 总 体 处 于 0.01 ~ 1 hm 之 
间 ,数量 约 占 斑 块 总 量 的 70% (图 3b)。 面积 低 于 
0.01 hm” 和 超过 10 hm 的 斑 块 较 少 ,数量 占 比 均 低 
F 10% (图 3b)。 城 市 建设 用 地 主要 分 布 在 正 西 和 
东北 方向 ,面积 均 超过 了 800 hm ;而 在 正 北 与 正 南 
方向 较 少 , 均 低 于 300 hm’ (图 4b)。 距 离 城市 中 心 
2 ~4 km 的 圈 层 是 城市 建设 用 地 分 布 最 集中 的 区 


(a) 城市 建设 用 地 提取 结果 


250 (b) 城市 建设 用 地 提取 结果 的 斑 块 分 布 特征 

< 196 
~ 200 180 
ni 
= 150 
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g 50 
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图 3 基于 全 卷 积 网 络 的 城市 建设 用 地 提取 结果 


Fig.3 Urban construction land extraction result based on fully convolutional network 
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(a) 方位 分 区 和 图 层 划 分 示意 图 
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4 2018 年 张掖 市 甘 州 区 建设 用 地 空间 分 布 特征 


Fig.4 Spatial patterns of urban construction land in Ganzhou district, Zhangye in 2018 


5 基于 KTH-Pavia 方法 的 城市 建设 用 地 提取 结果 


Fig.5 Urban construction land extraction result based on KTH-Pavia urban extractor 


域 ,超过 650 hm’; 而 距离 9 km 以 外 的 圈 层 城市 建设 
用 地 的 面积 很 少 ( 图 4c)。 


5 讨论 


5.1 新 方法 的 优势 和 特点 

在 现 有 方法 中 ,GAMBA 等 “ ”发 展 的 KTH- 
Pavia 方法 是 基于 Sentinel-1A 数据 提取 城市 建设 用 
地 信息 的 最 常用 方法 。 为 比较 全 卷 积 网 络 与 现 有 方 
法 的 提取 效果 ,首先 基于 Sentinel-1A 数据 ,使 用 
KTH-Pavia 方法 提取 出 了 研究 区 内 的 城市 建设 用 地 
(图 5) ,然后 在 提取 结果 中 对 城市 建设 用 地 与 非 城 
市 建设 用 地 像 元 进行 分 层 随 机 抽样 ,两 类 像 元 各 抽 
取 了 100 个 验证 点 进行 精度 评价 ,最 后 通过 构建 混 
淆 矩阵 后 计算 了 精度 评价 指标 ( 表 2) 。 

利用 全 卷 积 网 络 提取 的 城市 建设 用 地 结果 在 空 
间 形 态 上 更 接近 真实 地 物 。 提 取 结 果 形 态 完整 ,可 
以 清晰 ,准确 地 描绘 城市 建设 用 地 外 部 轮廓 ,能 够 很 
好 地 反映 城市 建设 用 地 的 整体 形态 (图 3a, 图 6a5 ) 。 


表 2 基于 KTH-Pavia 方法 的 城市 建设 用 地 
提取 结果 的 精度 评价 
Tab.2 Accuracy assessment of urban construction land 


extraction result based on KTH-Pavia urban extractor 


参考 图 像 


城市 建设 ” 非 城市 建设 总 和 
地 ( 像 元 ) 用 地 ( 像 元 ) ( 像 元 ) 


分 类 图 像 城市 建设 用 地 ( 像 元 ) 93 7 100 
非 城市 建设 用 地 ( 像 元 ) 31 69 100 
总 和 ( 像 元 ) 124 76 200 
总 体 精度 =81.00% Kappa =0.62 
制图 精度 ( 漏 分 误差 ) ”用 户 精 度 ( 错 分 误差 ) 
城市 建设 用 地 75.00% (25.00% ) 93.00% (7.00% ) 
非 城市 建设 用 地 90.79% (9.21% ) 69.00% (31.00% ) 


城市 内 部 建成 区 集中 连 片 的 分 布 形态 得 到 了 清晰 反 
映 ,提取 结果 空间 结构 紧凑 完整 (图 6b5 ) 。 利 用 
KTH-Pavia 方法 得 到 的 提取 结果 较为 零散 ,城市 建 
设 用 地 轮廓 不 清晰 ,该 问题 在 西部 区 域 尤为 明显 ,大 
量 城 市 建设 用 地 未 被 有 效 提取 (图 5, 图 6a7) 。 且 
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GS 泼 分 的 城市 
[|] 正确 分 类 的 非 城市 C EMAAR 
B 正确 分 类 的 城市 国 汪 正确 分 类 的 城市 


[| 错 分 的 城市 


|__| 正确 分 类 的 非 城市 [| 正确 分 类 的 非 城市 
GE 正确 分 类 的 城市 国 正确 分 类 的 城市 


上 错 分 的 城市 Maa 


E wit 


注 :范围 为 图 1 中 A B 两 个 缩放 窗口 。al 和 bl 为 真 彩 色 合 成 影像 ;2 和 b2 为 合成 孔径 雷达 VH 极 化 数据 ;a3 和 b3 为 合成 孔径 雷达 VV 极 化 


数据 ;a4 和 b4 为 目 视 解 译 获 取 的 真 值 ;a5 和 b5 为 全 卷 积 网 络 提取 结果 ,a6 和 b6 是 全 卷 积 网 络 提取 结果 与 真 值 的 对 比 ;a7 和 b7 为 KTH-Pavia 
提取 结果 ,a8 和 b8 是 KTH-Pavia 提取 结果 与 真 值 的 对 比 


图 6 结果 对 比 


Fig.6 Results comparison 


tió WN 


城市 建设 用 地 内 部 空间 结构 松散 , 漏 分 问题 严重 ,不 
能 有 效 反映 地 物 的 真实 空间 格局 (图 6b7) 。 利 用 全 
卷 积 网 络 提取 的 城市 建设 用 地 分 类 精度 更 高 ( 表 1， 
表 2) 。 利 用 全 卷 积 网 络 得 到 的 城市 建设 用 地 提取 
结果 Kappa 系数 达到 0. 85 ,利用 KTH-Pavia 方法 提 
取 的 城市 建设 用 地 Kappa 系数 为 0. 62 ,前 者 比 后 者 
高 出 约 37% 。 同 时 ,全 卷 积 网 络 提取 结果 的 总 体 精 
度 达 到 92.50% , 比 KTH-Pavia 方法 提取 结果 的 总 体 
精度 (81.00% ) 高 出 11% 以 上 。 全 卷 积 网 络 提取 结 
果 的 制图 精度 达到 89.72% ,而 KTH-Pavia 方法 仅 为 
75.00% ,前 者 比 后 者 高 出 近 15% 。 

基于 全 卷 积 网 络 的 城市 建设 用 地 提取 结果 拥有 
更 高 精度 的 主要 原因 ,是 全 卷 积 网 络 在 融合 和 提取 
Sentinel-1A 数据 多 极 化 类 型 与 多 尺度 纹理 信息 方面 
有 着 明显 优势 。 

首先 ,全 卷 积 网 络 支 持 不 同 层 数 和 大 小 的 数据 
输入 ,可 在 提取 过 程 中 充分 融合 不 同 极 化 方式 数据 
所 体现 的 地 物 散 射 特征 ,因此 提取 结果 更 加 接近 实 
际 地 物 ( 图 2b, 图 6al .6a4 .6a5 .6a6)。 而 KTH-Pavia 
在 图 像 信 息 提 取 过 程 中 无 法 对 数据 的 不 同 极 化 方式 
的 进行 融合 ,导致 局 部 提取 结果 仅 反 映 单 一 极 化 方 
式 的 纹理 特征 ,从 而 无 法 较 好 地 抑制 漏 分 问题 (图 
6a2 6a3 .6a7 .6a8 ) 。 

其 次 ,全 卷 积 网 络 包含 多 层 卷 积 结构 ,能 够 在 不 
同 尺 度 上 对 纹理 信息 进行 提取 ,可 充分 发 挥 Senti- 
nel-1A 数据 纹理 信息 丰富 的 优势 (图 2b)。 受 散射 
体 结构 影响 ,雷达 数据 会 在 局 部 出 现 异常 值 (如 建 
筑 物 的 “ 角 反 射 器 ”效应 会 产生 极端 高 值 ) ,而 全 卷 
积 网 络 对 多 尺度 信息 的 充分 融合 减弱 了 异常 值 引发 
的 负面 影响 ,从 而 有 效 抑 制 了 漏 分 问题 (图 6)。 卷 
积 核 的 作用 等 效 于 “激活 ”与 其 结构 相似 的 图 像 特 
征 , 多 尺度 卷 积 核 的 使 用 有 效 提 高 了 边缘 信息 与 纹 
理 特征 对 分 类 结果 的 贡献 中 。 经 过 充分 训练 后 ,全 
卷 积 网 络 可 以 有 效 表征 城市 地 物 不 同 尺 度 的 纹理 特 
征 ,充分 抑制 局 部 异常 值 影响 。 相 比 之 下 , KTH- 
Pavia 方法 则 依赖 灰 度 共生 和 矩阵 提取 单一 尺度 纹理 
言 息 ,导致 提取 结果 趋 于 破碎 并 出 现 明显 的 漏 分 问 
题 (图 6) 。 
5.2 展望 

本 研究 发 展 的 全 卷 积 网 络 不 仅 能 够 充分 发 挥 
Sentinel-1 A 数据 在 城市 建设 用 地 提取 方面 的 应 用 潜 
力 , 并 且 具 备 深度 学 习 方 法 的 迁移 学 习 能 力 , 可 以 应 
用 于 不 同 传感器 和 不 同 分 辩 率 的 数据 ,在 基于 诸如 


TerraSAR-X \Radarsat-2 等 不 同 SAR 数据 源 的 城市 
建设 用 地 提取 方面 具有 很 好 的 应 用 潜力 。 但 目前 的 
网 络 结构 仍 略 显 简单 , 且 所 使 用 的 训练 样本 数量 有 
限 ,使 得 模型 在 应 用 于 不 同 地 区 的 城市 建设 用 地 提 
取 时 可 能 尚 不 具有 理想 的 谤 化 能 力 。 

未 来 应 根据 不 同 地 区 的 Sentinel-1A 数据 特点 
对 全 卷 积 网 络 结构 进行 适当 微调 ,输入 不 同 代 表 性 
区 域 的 样本 对 网 络 进行 训练 ,以 提高 该 网 络 的 泛 化 
能 力 , 籍 此 则 可 利用 调整 后 的 网 络 在 更 大 尺度 上 实 
现 城市 建设 用 地 的 快速 高效 提取 。 


6 结论 


本 文通 过 构建 并 训练 全 卷 积 网 络 发 展 了 一 种 基 
于 Sentinel-1A 合成 孔径 雷达 数据 监测 城市 建设 用 
地 的 新 方法 ,使 用 该 方法 提取 了 张掖 市 甘 州 区 2018 
年 的 城市 建设 用 地 。 提 取 结 果 总 体 精 度 为 92. 50% ， 
Kappa 系数 为 0.85。 与 现 有 方法 中 应 用 最 广泛 的 
KTH-Pavia 方法 相 比 ,Kappa 系数 提高 了 37. 10% , 总 
体 精 度 提 高 了 11.50% 。 方 法 的 优势 主要 在 于 使 用 
全 卷 积 网 络 充分 融合 Sentinel-1A 数据 的 不 同 极 化 
方式 和 多 尺度 特征 ,从 而 提高 分 类 精度 。 

提取 结果 表明 研究 区 2018 年 城市 建设 用 地 面积 
为 4 006 hm’ , 斑 块 面积 集中 于 0.01 ~1 hm? 之 间 。 城 
市 主要 分 布 在 正 西 和 东北 方向 ,距离 城市 中 心 2~4 
km 的 圈 层 是 城市 建设 用 地 分 布 最 集中 的 区 域 。 

本 文 所 发 展 的 方法 不 仅 适 用 于 Sentinel-1A 数 
据 ,对 于 不 同 的 SAR 数据 源 均 具有 良好 的 应 用 洪 
力 。 此 外 ,通过 补充 训练 数据 和 微调 参数 , 即 可 迁移 
到 不 同 研究 区 使 用 ,拥有 较 强 的 可 迁移 性 。 
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A new method for detecting urban construction land based on Sentinel-1A 
synthetic aperture radar data and fully convolutional network 


LU Wen-lu'”, LIU Zhi-feng'*, HE Chun-yang'*, XIA Pei!” 
(1 Center for Human-Environment System Sustainability (CHESS) ,State Key Laboratory of Earth Surface Processes and 
Resource Ecology ( ESPRE) , Faculty of Geographical Science , Beijing Normal University , Beijing 100875 , China; 
2 School of Natural Resources , Faculty of Geographical Science , Beijing Normal University , Beijing 100875 , China ) 


Abstract: Timely and accurate extraction of urban construction land is essential for assessing the impacts of urban 
expansion on the environment. Optical remote sensing is susceptible to weather conditions and contains insufficient in- 
formation on the texture of urban features,and is thus not conducive to urban construction land extraction. Synthetic 
aperture radar (SAR) is unique in its potential to be used in all weather conditions and at any time. Furthermore , its 
backscattering and polarimetric information is sensitive to the dielectric and geometric properties of the urban land 
surface , providing a new means for rapid and accurate extraction of urban construction land. Among existing SAR da- 
ta, the Sentinel-1A SAR data is widely used due to its free access. Nonetheless , existing methods are poorly adapted to 
fully utilize the Sentinel-1 A SAR data , thereby limiting its application. The fully convolutional network is a deep learn- 
ing network developed on the basis of the convolutional neural network ,which adopts “pixel to pixel” image recogni- 
tion and has become an effective method for extracting urban construction land. This paper aims to develop a method 
for detecting urban construction land based on Sentinel-1A SAR data ,and thus a fully convolutional network. Firstly , 
a fully convolutional network consisting of five Inception modules was developed with reference to GoogleNet. Each In- 
ception module includes three convolutional layers with convolutional kernel sizes of 1 x 1,3 x3,and 5 x5 respective- 
ly ,three corresponding activation layers ,and a concatenation layer. The Ganzhou District in Zhangye City , Gansu Prov- 
ince ,China,where urban construction land and bare land show similar spectrum features , was used as a case study ar- 
ea to verify our fully convolutional network. The extraction results show that the area of construction land in the Gan- 
zhou District was 4,006 hectares in 2018, accounting for about 20% of the total area. The patch area of urban con- 
struction land was generally found to be between 0.01 and 1 hectare. Urban construction land was mainly distributed 
in the west and northeast ,and the circle of 2 -4 km from the city center exhibited the most concentrated urban con- 
struction land. Our accuracy assessment shows an overall accuracy of 92.50% and a Kappa coefficient of 0. 85. By 
comparing our results to extraction results using the “KTH-Pavia” method (the most widely used method for extrac- 
ting urban construction land from SAR data) ,it was found that extraction results based on the fully convolutional net- 
work are closer to the real urban construction land in spatial patterns. Furthermore ,the overall accuracy and Kappa co- 
efficient were ,respectively ,11% and 37% higher than the “KTH-Pavia” method. The principal reason for this higher 
accuracy is that the fully convolution network can better integrate multi-polarization and multi-scale texture informa- 
tion from SAR data. Furthermore , the fully convolutional network contains multiple convolutional structures and sup- 
ports multi-source data inputs. The method developed in this study has both greater accuracy than existing methods 
and is applicable to urban construction land extraction based on different SAR data in different regions. It therefore 
has potential for widespread application. 

Key words: Sentinel-1A; Sentinel data; synthetic aperture radar; urban construction land; fully convolutional 


network; deep learning 


